Entanglement Control of Azobenzene by Photoisomerization in 

NMR Quantum Computer 

Taiga Yasuda^, Masahito Tada-Umezaki^, Mikio Nakahara^'^ Tomonari Wakabayashi^'^ 

^Department of Physics, Kinki University, 
3-4- i Kowakae, Higashi- Osaka, 577-8502, Japan 
'^Research Center for Quantum Computing, 
Interdisciplinary Craduate School of Science and Engineering, 
Kinki University, 3-4-1 Kowakae, Higashi-Osaka 577-8502, Japan 
^Department of Chemistry, Kinki University, 
3-4-1 Kowakae, Higashi-Osaka, 511-8502, Japan 

Abstract 

Entanglement control of qubits in a photoisomerizing molecule is studied in the context of an 
NMR quantum computer by taking azobenzene as an example. Azobenzene has two different 
isomers, trans-azobenzene (TAB) and cis-azobenzene (CAB), which can be interconverted by pho- 
toisomerization. Changing molecular structure leads to change in the spin-spin coupling constant, 
and hence leads to change in entangling operation time. We first obtain stable structures of TAB 
and CAB by ah initio calculation. Then, we calculate the NMR spectra of these isomers and verify 
that they reproduce the chemical shift obtained experimentally with a good precision. Our result 
indicates that the coupling strength between a ^^N and a ^"^C nuclei in the molecule changes by a 
large amount under photoisomerization. 



1 



INTRODUCTION 



Quantum computer utilizes quantum mechanical phenomena, such as superposition 
and entanglement, as computational resources to outperform current digital computers. 
Nonetheless, a working quantum computer is yet to be realized. Feynman proposed the 
possibility of a quantum computer for the first time [l| and Deutsch triggered quantum 



computer research by introducing a quantum Turing machine [2j. Shor's algorithm |3[ and 
Grover's algorithm |^ demonstrated that a quantum computer can process some practi- 
cal problems more efficiently than a classical computer by making use of an exponentially 
massive parallel processing, known as the quantum parallelism. 

There are several physical realizations of a quantum computer proposed to date 
trapped ions neutral atoms in optical lattices js], an NMR quantum computer j^l, and 
superconducting qubits [lO|, among others. Nuclear spins in a molecule work as qubits, 
that are controlled by nuclear magnetic resonance (NMR) spectrometer in an NMR quan- 
tum computer ll|. An NMR quantum computer may have several qubits if an appropriate 
molecules are employed. ^^C-labelled chloroform 12| works as a two-qubit quantum com- 
puter, for example. 

Controlling the coupling between nuclear qubits in a molecule by structural change has 



been already attempted 



131 ]. In this paper, we analyze molecules whose inter-qubit coupling 



can be controllable by photoisomerization. Azo compounds are molecules whose two func- 
tional groups (aryl or alkyl) R and R' are connected by an azo group (N=N) as R-N=N-R'. 
They change their molecular structure between irans- azobenzene (TAB) and czs-azobenzene 
(CAB) under photoisomerization. Application of the photoisomerization of azobenzene to 



the switching device has been studied in 



14|-|l7|. In fact, CAB is unstable and it gradually 



changes to TAB, although the relaxation time of this process is long enough to control qubits 



many times [18|. One of the advantages of using azobenzene for quantum computation is 
that it is possible to control interaction between qubits by changing molecular structure by 
UV light irradiation. It is the purpose of this paper to study the interaction between nuclei 
in azobenzene and to propose the possibility of using azobenzene or similar molecules as a 
multi-qubit NMR quantum computer with a tunable couphng. 
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THEORETICAL BACKGROUND 



Let us consider a molecule with two spin-1/2 nuclei in a strong magnetic field Bo along 
the z-axis. We assume that the nuclear species of the spins are different. The Hamiltonian 
of the molecule in the laboratory frame is 

'H = -UJo,lIz® I -(^0,2! h+'Hint, (1) 

where /fc = (7^/2 is the kth component of the Pauli spin matrix divided by 2 and Wq.j is the 
Zeeman energy of the ith spin. The interaction Hamiltonian between the spins is given by 

^int= Yl Jjklj^h. (2) 

j,k=x,y,z 

Here Jjk are called the spin-spin coupling constants (SSCC). The coupling constants Jju 
take the isotropic Heisenberg form Jjk = J6jk for a molecule in a liquid state at room 
temperature, which we assume to be the case throughout this paper. Then "Hint simplifies 
as 

Hint = J Yl ^k®h. (3) 

k=x,y,z 

Let us fix the basis vectors as 



|0)= , |1)= , (4) 




where |0) (|1)) corresponds to the spin-up (spin-down) eigenstate of a^- We take the order 
of the two-qubit binary basis vectors as {|00), |01), |10), |11)} as usual. We transform the 
Hamiltonian in the laboratory frame ([1]) to that in the rotating frame of each spin by 
introducing the unitary transformation 

U{t) = e-'"""''^^' ® e-i^o,2ht_ (5^) 
The Hamiltonian in the rotating frame is 

-Hrot = unu^ - iu^u^ 

at 

= JIz®Iz. (6) 
where rapidly oscillating matrix elements have been dropped. 



3 



The time-evolution operator in this frame is 



f/rot(t) 



( ^-iJt/A 




V 





e*-^*/4 















-iJt/A 



\ 



(7) 



where we have taken the natural unit, in which h = 1. Suppose the initial state is a product 
state I + 0) = :^(|00) + |10)), for example, where |±) are the eigenvectors of ax with the 
eigenvalues ±1. This state may be generated by applying the Hadamard gate on the first 
qubit of a state |00), for example. The action of Urotit) on | + 0) yields the state 



^rot(t)|+0) 



^(e-'-'*/>0) + e^''*/^|10)) 
v2 



cos 



Jt 



0) 



I sm 



Jt 



0). 



(8) 



This shows that the spins are in a maximally entanged state (MES) at t = r such that 
I cos(Jr/4)| = I sin(Jr/4)|, that is. 



TT 



T=—(2n + l] 



[n 



0,1,2,...). 



(9) 



METHOD 



Stable structures and vibrational properties of azobenzene were already reported in |19| . 
in which GaussianOS with density functional theory (DFT) and M0ller-Plesset (MP) methods 
were employed. We used DFT method in the present work since MP method fails to calculate 
SSCC First, we evaluated optimized structures of azobenzene with DFT calculation: B3LYP 
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2lj, B3PW91 |2l|,|22|, PW91PW91 |22| and PBEPBE |23| with basis set: 6-31+G(d). 



Next, we calculated the NMR spectrum and SSCC of azobenzene using Gauge- 
Independent Atomic Orbital (GIAO) method 24N27l| and Integral-Equation-Formalism Po- 



larizable Continuum Model (lEFPCM) method [28|, |29|. GIAO method was employed to 
calculate the NMR spectrum, while lEFPCM method was used to take account of the 
solvent effect of chloroform in our calculation. Chemical shifts of carbons are relative to 
tetramethylsilane (TMS) and those of nitrogens are relative to ammonia (NH3). 
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RESULTS 



We calculated stable structures of azobenzene with GaussianOS with B3LYP/6-31+G(d), 
B3PW91/6-31+G(d), PW91PW91/6-31+G(d) and PBEPBE/6-31+G(d). The results ob- 
tained reproduce the previous ones [19| with a good precision (data not shown). Figured] 
shows the schematic structures of TAB and CAB. We studied entanglement between a ni- 
trogen nucleus and a carbon nucleus, which are denoted as nuclei 1 and 7 in Fig. [1] (a), 
and 1' and 7' in Fig. [1] (b), respectively. It should be noted that the N=N-C angle changes 
drastically by photoisomerization, under which the molecule transforms between TAB and 
CAB. We note from the NMR spectra of two nitrogen nuclei that they are very weakly 
coupled, that is J ^ for this pair, and hence it is takes extremely long time to entangle 
these nuclei. 

Table I shows the chemical shifts of the nitrogen and the carbon nuclei and SSCC between 
them we have obtained. They are computed with DFT functionals with B3LYP/6-31+G(d), 
B3PW91/6-31+G(d), PW91PW91/6-31+G(d) and PBEPBE/6-31+G(d). The times to at- 
tain the maximal entanglement were calculated with Eq. by setting n = 0. 

Experimentally measured values of the chemical shifts of the carbon and the nitrogen 
nuclei are also given in Table I 30|. The calculated chemical shifts of carbon nuclei reproduce 
the experimental results fairly accurately, while those of nitrogen nuclei do not agree with 
the experimental results with a good precision. This difference might be attributed to the 
difference in experimental conditions. The observed chemical shifts of carbon nuclei were 
obtained for azobenzene in liquid chloroform solvent, while those of nitrogen nuclei were 
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obtained for azobenzene in a polycrystalline state 

Experimental values of SSCC between the nitrogen and the carbon nuclei in azobenzene 
do not exist to our knowledge. In contrast, the measured data of SSCC between carbon 



nuclei are available 



32| and they are in good agreement with our computational results 



(data not shown). We expect, from these evidences, that the computed SSCC between the 
nitrogen and the carbon nuclei are reliable. Time to produce MES in CAB is longer than 
that of TAB in all computational results. The ratio of times required to attain MES in CAB 
to that of TAB depends on the scheme; it is 4.2, 3.6, 2.3 and 2.4 with B3LYP, B3PW91, 
PW91PW91 and PBEPBE, respectively. 

In this work, we have employed azobenzene to propose the possible application of pho- 
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toisomerizing molecules to an NMR quantum computer. The results showed that the time 
required to attain MES changes by a large amount between the two isomers, CAB and TAB. 
This suggests that photoisomerization can be used to control SSCC, and hence entanglement 
in azobenzene, which is potentially a powerful tool in NMR quantum computing. 

CONCLUSION AND DISCUSSION 

We calculated the time required to attain maximal entanglement between the carbon 
and nitrogen nuclei in azobenzene with ab initio methods. As a result, the time to attain 
maximal entanglement in czs-azobenzene is approximately four times longer than that in 
irans-azobenzene. This reflects the fact that SSCC of czs-azobenzene is four times weaker 
than that of irans-azobenzene. We conclude that azobenzene can be employed as a molecule 
with a tunable coupling in NMR quantum computing. 

Search for photoisomerizing molecules with a more drastic change in SSCC under pho- 
toisomerization is in progress and will be reported elsewhere. 
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(a) 

7' 



N = N 




3' 4' 



(b) 

FIG. 1: Schematic structures of trans-azohenzene (a) and cis- azobenzene (b). 



TABLE I: Chemical shift, spin-spin couphng constant and time to attain maximahy entangled state 
of irans-azobenzene (upper table) and cis-azobenzene (lower table). '^Chemical shift of nitrogen 
nucleus (7 and 7') and '^chemical shift of carbon nucleus (1 and 1'). Experimental values are from 



=Ref. (3l|], '^Ref. [30|1, and ''Ref. [3^. The value of J(12, 16) is an average of calculated J(12) and 
J(16) since the experimentally available value is an average of J(12) and J(16). This also applies 
to J(1'2',1'6'). 
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